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RAND (Computation of Random Fractions) 


DECUS Program Library Write-up 


DECUS No. 8-319 


ABSTRACT 

I h X° 9ra T T COmpufes a des 'gnated number of random fractions in the range of 0.0000 to 
U.V999. The computed fractions represent a sample from a uniform distribution. 

TAPES REQUIRED 


1 ‘ . f orf n of program tape - The program is written in the PDP-8 FORTRAN-D lanauaae 
and is in the source language. ' 

OPERATING INSTRUCTIONS 

Normal for the disk operating system. On continuing after the teletype has printed READY 
he program pauses for the entry of five integers. The first four integers must consist of two ' 
digits each, and are used to specify the starting point for the series of pseudo-random 
numbers computed. The fifth integer specifies the number of random fractions required. 


53 95 27 04 40 

from n d r/ he ke ? b r ( itl™? ef 7 pe WOuld request the prinHn 9 of 40 ™dom fractions 
from a starting point of 53952704. (Note that a given starting point will always generate 

he same sequence of random fractions, a feature which is frequently of value in the use of 

random fractions in simulation.) A result tape can be obtained by switching on the low- 

speed punch before keying "return" after entering the last of the integers. 

OUTPUT 

The program prints random fractions in the range of 0.0000 to 0.9999. The fractions are 

printed four to a line, the total number of fractions being specified by the fifth of the control 
integers. 


STORAGE AND LIMITATIONS 
Normal for FORTRAN-D 
METHOD 


The program uses Wingersky’s FORTRAN function subprogram which implements the power 

u!e of e i„ m re .tf"!™?" 9 pseudo-random numbers. This subprogram incorporates an ingenious 

method tT Pr ° 9ram ,he modul ° ' multiplication necessary for the power residue 

^cedures! '"wby 1968 "* L ° h " eS ° nd C 00 '^' 5 "'"Action *> Statistical 
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C PROGRAM TO COMPUTE RANDOM FRACTIONS 

L=-l 

ACCEPT 3, K1, K2, K3, K4, N 
3 FORMAT (1/1/1/1/0 

DO 10 1=1, N 
M1=11*K4 
M2=l 1 * K3 
M3=l 1*K2+K4 
M4=11*K1+K3 

J=M1/100 

K4=M1-100*J 

M2=M2+J 

J=M2/100 ^ 

K3=M2-100*J 

M3=M3+J 

J=M3/100 

K2=M3-100*J 

M4=M4+J 

J=M4/100 

K1=M4-100*J 

X1=K1 

X2-K2 

R=X1 *1.0E-2+X2*l .0E-4 
L=L+1 

IF (L—4) 1,2,2 
2 1=0 

TYPE 101, R 
GO TO 10 

1 TYPE 102, R % 

10 CONTINUE 

101 FORMAT (/,E) 

102 FORMAT (E) 

STOP 

END 
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